Dipole formation at metal/PTCDA interfaces: role of the Charge Neutrality Level 
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The formation of a metal/PTCDA (3, 4, 9, 10- perylenetetracarboxylic dianhydride) interface 
barrier is analyzed using weak chemisorption theory. The electronic structure of the uncoupled 
PTCDA molecule and of the metal surface is calculated. Then, the induced density of interface 
states is obtained as a function of these two electronic structures and the interaction between both 
systems. This induced density of states is found to be large enough (even if the metal/PTCDA 
interaction is weak) for the definition of a Charge Neutrality Level for PTCDA, located 2.45 eV above 
the highest occupied molecular orbital. We conclude that the metal/PTCDA interface molecular 
level alignment is due to the electrostatic dipole created by the charge transfer between the two 
solids. 

<N 

Electronic materials made of molecular films are a fast developing field, with many potential applications in organic- 
based devices. Designing new organic-based materials requires a detailed understanding of the different processes 
occurring in these devices. In particular, metal/organic and semiconductor/organic interface barriers play a decisive 
i— h ■ role 0,0. However, the formation of barriers is not yet well understood. 

In the Schottky-Mott model of metal/organic interfaces, it is assumed that no interface dipole is formed at the 
junction, and that the position of molecular levels with respect to the metal Fermi level is defined by vacuum level 
alignment. This situation was disproved by Narioka et al. [j| who, using ultra-violet photoemission spectroscopy 
(UPS), found large interface dipoles (~ 0.5 — l.OeV) at several metal/organic interfaces. Independent data by Hill et 
al. |4j confirmed this conclusion. Various mechanisms are believed to operate simultaneously at these interfaces, and 
I \ several models have been advanced [lj|2| ■ Metal- molecule chemical reaction has been seen to create interface gap states 
that pin the Fermi level Q , a situation that is analogous to that described by the Unified Defect Model proposed for 
inorganic semiconductor/metal interfaces [||. Compression of the metal surface electronic tail by adsorbed molecules, 
leading to vacuum level interface shift (the "pillow" effect), has also been proposed as a general metal/organic interface 
mechanism @,HI3- 

In this letter, we explore the first application to a metal/organic interface of the Induced Density of Interface 
(or virtual) States (IDIS) Model \u%. We study a metal/PTCDA (3, 4, 9, 10- perylenetetracarboxylic dianhydride) 
\ interface and analyze how the chemical interaction between the organic molecule and the metal creates an IDIS in the 
^v^j ' organic energy gap. Our calculations show that, although the chemical interaction is weak, the IDIS is large enough 
\ that a Charge Neutrality Level (CNL) of the organic molecule can be defined. Our results show that the interface 
. Fermi level Ep is pinned at the CNL, a situation similar to that described for the formation of Schottky barriers at 
' conventional semiconductor/metal junctions. 

In this theoretical analysis, we study the metal/PTCDA interaction in three steps. First, the molecular orbitals of 
the organic molecule are calculated using a Density Functional Theory (DFT) local-orbital method [ill . Care is 
taken to identify the appropriate wavefunctions associated with the Affinity and Ionization levels defining the molecular 
transport gap. Then, we calculate the metal/molecule interactions and obtain the electronic induced density of states 
in the molecular energy gap. In this approach, we neglect the molecule-molecule interaction, which introduces only 
a small broadening of the molecular levels (~ 0.2 eV), and does not create any electron density of states in the 
molecular energy gap. Finally, we define the CNL associated with these induced density of states and show that the 
J-^ ' organic/metal interaction creates density of states large enough to pin the Fermi level close to the CNL. Pinning can 
^ , be characterized by the parameter 

• i-H , 

X- 

£ : S = §^, (1) 

(where 4>m is the metal workfunction) , equivalent to the interface slope parameter in the theory of 
metal/semiconductor junctions 13] and which we find here to be ~ 0.13. 

In our DFT study of the PTCDA molecule [l2^. we use an optimized minimal local orbital basis, and introduce 
many-body effects by means of a local orbital exchange and correlation potential 11]. In our actual calculations, we 
combined a local-orbital DFT code pjj with this exchange and correlation potential. This approach is equivalent to 
conventional DFT but is especially advantageous in cases where correlation is important 0] , and for calculating the 
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FIG. 1: DFT spectrum (long bars: n states; short bars: a states) and DOS (introducing a 0.5 eV FWHM Gaussian broadening) 
for the PTCDA molecule. 



Affinity (A) and Ionization (/) levels of the molecule, as will be seen below in the discussion of Koopman's theorem 

Figure^shows the spectrum obtained from the DFT calculation; we also calculate a LDOS by introducing for each 
level a Gaussian broadening with a 0.5 eV FWHM 16]. This DFT calculation produces a molecular energy gap ~ 1.6 
eV, very small compared to the theoretical single-particle gap A — I (~ 5eV) of the isolated molecule [13, and even 
smaller than its optical gap, ~ 2.6eV 2]. This is due to the fact that the DFT eigenvalues are not directly related to 
the molecular levels. Note also that, in analyzing the metal/organic interface, one needs to consider the single-particle 
energy gap of the isolated molecule, A — I. At the interface, this energy gap, A — I, is reduced by solid-state effects 
0, EI EH j as wm be discussed later. 

We first calculate A and / using the definitions A = E[N+1] —E[N], and I = E[N] - E[N - 1], where E[Ni\ is the 
ground-state energy of the molecule with Ni electrons. Each case can be calculated using the DFT method described 
above; we find [lj / = —5.9 eV and A = —1.1 eV, yielding an energy gap of 4.8 eV, in reasonable agreement with 
the value of ~ 5 eV 0|. Interestingly, similar results can be obtained using a sort of Koopman's theorem which 
amounts to introducing an electron (hole) in the LUMO (HOMO) of the neutral molecule, and calculating the total 
energy of the system neglecting electron relaxation effects. In this approach, the Affinity level, for instance, is simply 
the LUMO level corrected by a quantity 8 A ^1 > which can be interpreted as a self-interaction ener gy a ssociated with 
the new electron charge at the LUMO wavefunction. Likewise, the ionization-level correction, SI 12], is a negative 
self-interaction energy related to the removal of one electron from the HOMO. 

This Koopman's theorem approach yields I = — 6.1eV and A = —1.2eV, in very good agreement with the values 
given above. This also shows that electron relaxation effects in PTCDA associated with introducing or extracting an 
electron from the molecule are small, and indicates that the molecular orbitals calculated within DFT for the HOMO 
and LUMO levels are very appropriate for describing the Ionization and Affinity wavefunctions of the molecule. Note 
that 7T and a orbitals give different values of 8 A and 81. In our calculations, we obtain the PTCDA energy levels for 
the molecule by modifying the DFT tt and a levels with the corresponding values of 8 A and 81 given by Koopman's 
theorem. 

Solid-state effects associated with long-range electronic polarization have been analyzed by Tsiper et al. |18| . 
Polarization of molecules of the solid screens the electric field created by the extra charge introduced in the system, 
either an electron or a hole, modifying the energy levels I and A and the corresponding single-particle energy gap (also 
called transport gap in the solid). This correction is important, 1 — 1.5eU, yet the ionization and affinity wavefunctions 
of the molecule are not expected to present important modifications. Other effects, such as lattice relaxation and 
vibronic coupling, only introduce corrections of the order of 0.2eV, in the transport gap [l6j . 

Thus, the electronic wavefunctions of PTCDA can be described using the results of a DFT approach, introducing 
a Koopman's correction for the tt and a molecular energy levels. Polarization effects are then included by fitting the 
ionization level and transport gap to their experimental values by a rigid shift in the occupied and empty electron 
states. We fit the energy gap, A — I, to 3.2 eV, which seems to include all the effects discussed above. 

We analyze the metal/PTCDA contact using chemisorption theory in the limit of weak interaction. The inset 
of Figure |21 shows the geometry of interest: a planar PTCDA molecule located at a distance d from the outermost 
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FIG. 2: C, O and H - Au interaction as a function of PTCDA-Au distance. 



metal layer, i.e. the Au (111) surface. We assume the PTCDA wavef unctions, Vi> to be described by the DFT 
method discussed above, and the metal wavefunctions by a DFT-LDA local-orbital code 0], which also yields the 
Green-function, G a p(E), where a and (5 refer to the local-orbital basis used in the DFT code. In the limit of weak 
metal-PTCDA interacion ^(|, the main effect of the metal on the PTCDA electronic structure is to introduce the 
following self-energy for the molecular levels: 



(2) 



where is the hopping interaction between the molecular orbital ipi and the metal eigenfunction, ip v , G V (E) being 
the corresponding Green-function. Equation [21 can be rewritten in a more appropriate way by using a local-orbital 
basis for the molecule and the metal. With the notation ipi = J2j Cij^ji equation takes the form 
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(3) 



where the self-energy of each molecular level is determined by: (a) the hopping interaction Tj a between local 
orbitals; (b) the coefficients Cjj of the molecular orbitals ipi; and (c) the Green-function of the metal, G a p(E). 

Figure[2]shows the hopping interaction between the Au 6s orbital and the different orbitals of PTCDA, C 2s and 2p, 
O 2s and 2p, and H Is (in the following, we neglect the 5d and 6p orbitals of Au, which is a reasonable approximation 
because we are interested in the induced density of states around the metal Fermi energy) . These hopping elements 
are calculated using DFT for each pair of molecule-substrate atoms. 

Equation can be simplified to a more convenient form by neglecting contributions having j ^ j'\ these terms 
represent interference effects between different local orbitals of PTCDA, and tend to cancel each other out. The 
contribution coming from the Au 6s orbitals yields 



Z t (E) =J2\^j\ 2 \ T 3fis\ 2 G 6s (E). 



(4) 



Each molecular level contributes to the LDOS with —— Im < ■ 

7T I ( 



{e~e )-t. (e) } ■ T ne mam effect of Sj is to broaden 
the molecular levels by Im(Tn{E)). Since we are studying the density of states around the molecular energy gap, we 
take Ti/2 ~ 'E(Ep) {Ep is the interface Fermi level), and obtain 
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FIG. 3: PTCDA LDOS for the PTCDA/Au(lll) interface (top view at inset) calculated for d = 2.8 A (dashed line), d = 3.0 A 
(continuous line), and d = 3.2 A (dotted line). Long (short) bars correspond to the 7r (a) states neglecting the metal-molecule 
interaction. 



for the contribution of each molecular level to the LDOS. 

As shown in Figure|3 Tjfi s and Ti depend on the PTCDA- Au distance, d. We do not know of any experimental value 
for this interfacial distance, nor are DFT calculations reliable in obtaining it, since the metal/PTCDA interaction has 
an important Van der Waals (or dispersion forces) contribution. Recently, Krause et al. |2£| measured the PTCDA- 
Ag(lll) distance to be ~ 2.85 A, with an uncertainty of ±0.2 A. We expect a similar (or slightly larger) distance, of 
- 3.0 A for the PTCDA- Au(lll) system. 

Figure shows the PTCDA LDOS (both spins included) around the organic energy gap for d = 2.8,3.0 and 3.2 
A. For these distances, the tt and a levels are broadened by TJ ~ 2.1, 1.5 and 1.0 eV and Tf ~ 1.05,0.75 and 0.5 
eV, respectively. The interactions Tj^s and the level broadening terms Tf and Tf increase with decreasing distance: 
for d — 2.8 A the values of T are typically twice the values for d — 3.2 A. Figure |31 also shows the position of the 
CNL. The CNL is calculated by integrating the induced LDOS and imposing charge neutrality conditions: the total 
number of electrons up to the CNL equals that of the isolated molecule. The CNL is practically the same for the 
three distances d — 2.8, 3.0 and 3.2 A, showing that it is a very robust quantity, practically independent of the details 
of the metal/PTCDA interaction. 

We calculated the interface slope parameter, S — gf^-- In our case fl3j ]. 



d(j) M l + 47re 2 L>(£ F )<5/A' 1 ' 

where D(Ep) is the PTCDA LDOS (Figure[3J) at the Fermi level (practically, the CNL), 8 is the distance between 
the charges induced in the metal and the organic molecule (between 2.8 — 3.2 A in Figure|3J) and A is the area associated 
with a PTCDA molecule (around 120 A 2 ). We hnd a DOS of 1.5, 1.2 and 0.9 x 10 14 eV^ 1 cm' 2 for d = 2.8, 3.0 and 3.2 
A, respectively. For these distances, S ~ 0.12,0.13 and 0.16. At variance with the CNL, S strongly depends on the 
metal-PTCDA distance. But even for d = 3.5 A, we obtain a DOS of 0.5 x 10 14 eV -1 cm -2 and S ~ 0.23, indicating 
that the IDIS is still significant. This small value of S (0.1 — 0.2) shows that the LDOS induced in the molecule is 
high enough (even if the metal/PTCDA interaction is small) to strongly pin the interface Fermi energy close to the 
CNL. 

Up to this point we have assumed the molecule to present no geometrical distortion when deposited flat on the 
metal. Recent DFT calculations for organic molecules adsorbed on metals [2lJ show non-negligible bond-length 
changes (e.g. the C-C bond in an octane molecule adsorbed on Cu(110) is shortened by ~ 0.04 A). However, given 
the linear geometry of the alkanes, their geometrical distortions are probably larger than the ones in PTCDA on 
Au(lll). Nevertheless, we explored the sensitivity of our results to geometrical distortions by changing some of the 
molecular bond lengths. We analyzed how the deformation of the molecular C-C and C-0 bonds modify the position 
of the CNL and the LDOS. Our calculations indicate that the CNL shifts by 0.1 eV for a change of 0.1 A of the 
carboxyl C-C bonds and by ~ 0.07 eV for a 0.1 A change in the C-0 bonds, while in all cases the LDOS around the 
CNL is practically unaltered. Thus, we conclude that the CNL undergoes only minor changes due to the molecular 
deformation, and that the LDOS is not significantly modified. The CNL shift due to molecular distortions can be 
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estimated by realizing that the weaker single C-C bonds are probably going to suffer the largest deformation. In 
the case of the alkanes, these distortions are ~ 0.05 A 21]. This suggests that the accuracy in our calculated CNL, 
neglecting molecular deformations, is better than ±0.1 eV. 

From our calculations and the previous discussion, we find that the CNL is located 2.45 ±0.1 eV above the center of 
the HOMO level of the molecule, in very good agreement with the experimental Ep pinning position 4]. The value of 
S ~ 0.13 is also in good agreement with S ~ inferred from experiments The important outcome of our analysis 
is that the calculated value of S ~ 0.13 is indicative of a significant DOS at the metal/PTCDA interface, in spite of 
the weak chemical interaction between the two materials (see Figure [2J. The implication is that the interface Fermi 
level is close to the CNL (in our calculations, Ep lies ~ 0.02 eV below the CNL). Thus, the formation of the interface 
barrier is related to the transfer of charge across the interface. The charge transfer is associated with the IDIS, and 
creates an electrostatic interface dipole, which tends to align the metal Fermi level and the PTCDA CNL. 

Recently, Bagus et al. [2^ have suggested that for molecules physisorbed on metal surfaces a significant interface 
dipole, ~ 0.3 eV, originates from exchange-like effects (this is the "pillow" effect mentioned above 0,11; @)- This 
effect can be easily incorporated in our discussion as a modification of the initial metal workfunction. As the final 
interface Fermi level Ep is given by: 

{Ep - CNL) = S(4> M - CNL), (7) 

where Ep, CNL and 4>m must have a common origin, a change of 0.3 eV in 4>m is reflected in a Fermi level shift of 
0.04 eV (S = 0.13). This suggests that for PTCDA, where S is small, the "pillow" effect is almost negligible with 
regard to the final Fermi level position. 

In conclusion, we have analyzed the metal/PTCDA interface barrier formation using weak chemisorption theory. 
The aim of this paper was to explore the importance of the IDIS and its contribution to the interface formation. Our 
results show that even for Au-PTCDA distances for which the interaction is weak (d ~ 3.0 — 3.5 A), the IDIS is high 
enough to create an interface dipole that tends to align the Fermi level and the PTCDA CNL. A measure of this drive 
to align the two levels is given by the slope parameter, S, which our calculations show to be ~ 0.13. 

Within this framework, we conclude that the main mechanism that leads to the formation of the Schottky barrier in 
metal/PTCDA interfaces is the charge transfer associated with the tunneling of metal electrons through the molecular 
energy gap. Our results show that the interface Fermi level should be close to the CNL of the molecule, which our 
calculations place 2.45 ± 0.1 eV above the HOMO level. 
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